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MAGNETIC RESONANCE IMAGING OF 
PROSTATE BRACHYTHERAPY SEEDS 

BACKGROUND OF THE INVENTION 

[0001 ] The field of the invention is nuclear magnetic resonance imaging (MRI) 

methods and systems. More particularly, the invention relates to the imaging of soft 
tissues and implanted devices such as brachytherapy seeds. 
[0002] When a substance such as human tissue is subjected to a uniform 

magnetic field (polarizing field B 0 ), the individual magnetic moments of the spins in 
the tissue attempt to align with this polarizing field, but precess about it in random 
order at their characteristic Larmor frequency. If the substance, or tissue, is subjected 
to a magnetic field (excitation field B,) which is in the x-y plane and which is near the 
Larmor frequency, the net aligned moment, M z , may be rotated, or "tipped", into the x- 
y plane to produce a net transverse magnetic moment M t . A signal is emitted by the 
excited spins after the excitation signal B, is terminated, this signal may be received 
and processed to form an image. 

[0003] When utilizing these signals to produce images, magnetic field 

gradients (G x G y and G 2 ) are employed. Typically, the region to be imaged is scanned 
by a sequence of measurement cycles under the direction of prescribed pulse sequences 
in which these gradients vary according to the particular localization method being 
used. The resulting set of received NMR signals are digitized and processed to 
reconstruct the image using one of many well known reconstruction techniques. 
[0004] An advantage of MR imaging over other imaging modalities such as x- 

ray CT is its ability to differentiate between different soft tissues. The amplitude of the 
received NMR signal during a measured cycle is determined by a number of factors 
related to the molecular bonding of the excited spins with surrounding atoms. These 
factors can be magnified, or "weighted" by the choice of imaging pulse sequence scan 
parameters to enhance the contrast between different tissue types. 
[0005] Unfortunately, MR imaging does not provide optimal contrast between 

certain implants and surrounding soft tissues. Seeds used in brachytherapy of the 
prostate, for example, are small cylindrical-shaped objects made of titanium which are 
implanted in a malignancy. Accurate placement of these seeds is critical for the proper 
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treatment of the disease and it is normal practice to image the prostate following 
implantation to confirm placement. Currently, x-ray CT is the preferred imaging 
modality for this purpose because it produces an image in which the seeds are highly 
contrasted with surrounding soft tissues. The limitation with x-ray CT, however, is 
that it does not contrast one soft tissue from another very well, making it difficult to 
discern the prostate tissue and the malignancy therein. 

SUMMARY OF THE INVENTION 

[0006] The present invention is a method for acquiring and reconstructing an 

MR image which provides the high contrast between different tissues normally 
associated with this imaging modality and high contrast between implanted objects and 
surrounding soft tissues. More particularly, a complex k-space data set is acquired 
with the MRI system; a magnitude image is reconstructed from the k-space data set; a 
phase image is reconstructed from the k-space data set; the location of the implanted 
objects is determined from the phase image; and the determined location of the 
implanted objects is indicated on the magnitude image. 

[0007] An obj ect of the invention is to produce an image which differentiates 

between soft tissues and differentiates between implanted objects and surrounding soft 
tissues. It has been discovered that a phase MR image provides sharp contrast between 
implants and surrounding soft tissue. When combined with a conventional magnitude 
image that differentiates between soft tissues, the phase image provides the additional 
contrast mechanism for locating the implants. 

[0008] The foregoing and other objects and advantages of the invention will 

appear from the following description. In the description, reference is made to the: 
accompanying drawings which form a part hereof, and in which there is shown by way 
of illustration a preferred embodiment of the invention. Such embodiment does not 
necessarily represent the full scope of the invention, however, and reference is made 
therefore to the claims and herein for interpreting the scope of the invention. 



BRIEF DESCRIPTION OF THE DRAWINGS 

[0009] Fig. 1 is a block diagram of an MRI system which employs the present 

invention; 
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[0010] Fig. 2 is an electrical block diagram of the transceiver which forms part 

of the MRI system of Fig. 1; 

[001 1 ] Fig. 3 is a graphic representation of a gradient-recalled echo pulse 

sequence used in the MRI system of Fig. 1 to practice the preferred embodiment of the 
invention; 

[0012] Fig. 4 is a graphic representation of a spin-echo pulse sequence used in 

the MRI system of Fig. 1 to practice the preferred embodiment of the invention; and 
[0013] Fig. 5 is a flow chart of the steps performed by the MRI system of Fig. 

1 to practice the preferred embodiment of the invention. 

DESCRIPTION OF THE PREFERRED EMBODIMENT 
[0014] Referring first to Fig. 1 , there is shown the major components of a 

preferred MRI system which incorporates the present invention. The operation of the 
system is controlled from an operator console 100 which includes a keyboard and 
control panel 102 and a display 104. The console 100 communicates through a link 
116 with a separate computer system 107 that enables an operator to control the 
production and display of images on the screen 104. The computer system 107 
includes a number of modules which communicate with each other through a 
backplane. These include an image processor module 106, a CPU module 108 and a 
memory module 113, known in the art as a frame buffer for storing image data arrays. 
The computer system 107 is linked to a disk storage 1 1 1 and a tape drive 1 12 for 
storage of image data and programs, and it communicates with a separate system 
control 122 through a high speed serial link 115. 

[0015] The system control 122 includes a set of modules connected together by 

a backplane. These include a CPU module 1 19 and a pulse generator module 121 
which connects to the operator console 100 through a serial link 125. It is through this 
link 125 that the system control 122 receives commands from the operator which 
indicate the scan sequence that is to be performed. The pulse generator module 121 
operates the system components to carry out the desired scan sequence. It produces 
data which indicates the timing, strength and shape of the RF pulses which are to be 
produced, and the timing of and length of the data acquisition window. The pulse 
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generator module 121 connects to a set of gradient amplifiers 127, to indicate the 
timing and shape of the gradient pulses to be produced during the scan. The pulse 
generator module 121 also receives patient data from a physiological acquisition 
controller 129 that receives signals from a number of different sensors connected to the 
patient, such as ECG signals from electrodes or respiratory signals from a bellows, 
And finally, the pulse generator module 121 connects to a scan room interface circuit 
133 which receives signals from various sensors associated with the condition of the 
patient and the magnet system. It is also through the scan room interface circuit 133 
that a patient positioning system 134 receives commands to move the patient to the 
desired position for the scan. 

[00 1 6] The gradient waveforms produced by the pulse generator module 1 2 1 

are applied to a gradient amplifier system 127 comprised of G x , G y and G z amplifiers. 
Each gradient amplifier excites a corresponding gradient coil in an assembly gener ally 
designated 139 to produce the magnetic field gradients used for position encoding 
acquired signals. The gradient coil assembly 139 forms part of a magnet assembly 141 
which includes a polarizing magnet 140 and a whole-body RF coil 152. A transceiver 
module 150 in the system control 122 produces pulses which are amplified by an RF 
amplifier 151 and coupled to the RF coil 152 by a transmit/receive switch 154. The 
resulting signals radiated by the excited nuclei in the patient may be sensed by the 
same RF coil 1 52 and coupled through the transmit/receive switch 1 54 to a 
preamplifier 153. The amplified NMR signals are demodulated, filtered, and digitized 
in the receiver section of the transceiver 150. The transmit/receive switch 154 is 
controlled by a signal from the pulse generator module 121 to electrically connect the 
RF amplifier 151 to the coil 152 during the transmit mode and to connect the 
preamplifier 153 during the receive mode. The transmit/receive switch 154 also 
enables a separate RF coil (for example, a head coil or surface coil) to be used in either 
the transmit or receive mode. 

[001 7] The NMR signals picked up by the RF coil 1 52 are digitized by the 

transceiver module 150 and transferred to a memory module 160 in the system control 
122. When the scan is completed and an entire array of data has been acquired in the 
memory module 160, an array processor 161 operates to Fourier transform the data 
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into an array of image data. This image data is conveyed through the serial link 1 15 to 
the computer system 1 07 where it is stored in the disk memory 111. In response to 
commands received from the operator console 100, this image data may be archived on 
the tape drive 1 12, or it may be further processed by the image processor 106 and 
conveyed to the operator console 100 and presented on the display 104. 
[00 1 8] Referring particularly to Figs. 1 and 2, the transceiver 1 50 produces the 

RF excitation field Bl through power amplifier 151 at a coil 152 A and receives the 
resulting signal induced in a coil 152B. As indicated above, the coils 152A and B may 
be separate as shown in Fig. 2, or they may be a single wholebody coil as shown in 
Fig. 1 . The base, or carrier, frequency of the RF excitation field is produced under 
control of a frequency synthesizer 200 which receives a set of digital signals (CF) from 
the CPU module 1 19 and pulse generator module 121. These digital signals indicate 
the frequency and phase of the RF carrier signal produced at an output 201 . The 
commanded RF carrier is applied to a modulator and up converter 202 where its 
amplitude is modulated in response to a signal R(t) also received from the pulse 
generator module 121. The signal R(t) defines the envelope of the RF excitation pulse 
to be produced and is produced in the module 121 by sequentially reading out a series 
of stored digital values. These stored digital values may, in turn, be changed from the 
operator console 100 to enable any desired RF pulse envelope to be produced. 
[00 1 9] The magnitude of the RF excitation pulse produced at output 205 is 

attenuated by an exciter attenuator circuit 206 which receives a digital command, TA, 
from the backplane 118. The attenuated RF excitation pulses are applied to the power 
amplifier 151 that drives the RF coil 152A. For a more detailed description of this 
portion of the transceiver 122, reference is made to U.S. Patent No. 4,952,877 which is 
incorporated herein by reference. 

[0020] Referring still to Fig. 1 and 2 the signal produced by the subject is 

picked up by the receiver coil 152B and applied through the preamplifier 153 to the 
input of a receiver attenuator 207. The receiver attenuator 207 further amplifies the 
signal by an amount determined by a digital attenuation signal (RA) received from the 
backplane 118. 
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[002 1 ] The received signal is at or around the Larmor frequency, and this high 

frequency signal is down converted in a two step process by a down converter 208 
which first mixes the NMR signal with the carrier signal on line 201 and then mixes 
the resulting difference signal with the 205 MHz reference signal on line 204. The; 
down converted NMR signal is applied to the input of an analog-to-digital (A/D) 
converter 209 which samples and digitizes the analog signal and applies it to a digital 
detector and signal processor 210 which produces 16-bit in-phase (I) values and 16-bit 
quadrature (Q) values corresponding to the received signal. The resulting stream of 
digitized I and Q values of the received signal are output through backplane 1 18 to the 
memory module 1 60 where they are employed to reconstruct an image. 
[0022] The 2.5 MHz reference signal as well as the 250 kHz sampling signal 

and the 5, 10 and 60 MHz reference signals are produced by a reference frequency 
generator 203 from a common 20 MHz master clock signal. For a more detailed 
description of the receiver, reference is made to U.S. Patent No. 4,992,736 which is 
incorporated herein by reference. 

[0023] To practice the preferred embodiment of the invention two different 

pulse sequences are performed by the MRI system of Fig. 1 . One of these is a 
gradient-recalled echo pulse sequence shown in Fig. 3 and the other is a spin echo 
pulse sequence shown in Fig. 4. 

[0024] The gradient-recalled pulse sequence in Fig. 3 performs a slice selection 

by applying a selective RF excitation pulse 220 in the presence of a z-axis slice select 
gradient waveform 222 to produce transverse magnetization. An NMR echo signal 
230 is produced and acquired in the presence of an x-axis readout gradient 228 which 
encodes for position along a readout axis. This echo signal 230 is referred to as a 
gradient-recalled NMR echo signal because the transverse magnetization is first 
dephased by a negative lobe at the beginning of the readout gradient waveform 222 
and then re focused by the positive readout gradient lobe. The NMR echo signal 230 is 
position encoded along a y-axis by a phase encoding gradient pulse 226. The phase of 
the spin magnetization is subsequently rewound by a corresponding rewinder gradient 
pulse 234 of the opposite polarity directed along the phase encoding axis. During a 
scan the pulse sequence is repeated with a different phase encoding which is 
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incremented through 256 discrete values (-128 to +128). The particular scan 
parameters used will depend on the particular clinical application, but when applied to 
imaging the prostate an echo time (TE) of as low as 4 ms, a transmit repeat (TR) time 
ranging from 300 to 400 ms, and a receiver bandwidth of 16 kHz is preferred. 
[0025] The spin-echo pulse sequence of Fig. 4 acquires an image of exactly the 

same size and location as the gradient-recalled echo pulse sequence of Fig. 3. It 
performs a slice selection by applying a 90° selective RF excitation pulse 300 in the 
presence of a z-axis gradient pulse 301 and its associated rephrasing pulse 302. After 
an interval TE/2, a 180° selective RF pulse 303 is applied in the presence of another z- 
axis gradient pulse 304 to refocus the transverse magnetization at the time TE, and 
produce an echo NMR echo signal 305. The NMR signal 305 is referred to as a "spin- 
echo" NMR signal because it is produced by transverse magnetization which is 
refocus ed by an RF refocusing pulse 303. 

[0026] To position encode the spin-echo signal 305, an x-axis read-out gradient 

pulse 306 is applied during the acquisition of the NMR signal 305. The read-out 
gradient frequency encodes the NMR signal 305 in the well known manner. In 
addition, the echo signal 305 is position encoded along the y-axis by a phase encoding 
gradient pulse 307. The phase encoding gradient pulse 307 has one strength during 
each pulse sequence, and it is incremented in steps through 256 discrete strengths (-128 
to +128) during the scan. As a result, each of the 256 NMR spin-echo signals 305 
acquired during the scan is uniquely phase encoded. 

[0027] As with the gradient-recalled echo pulse sequence, the particular scan 

parameters for this spin-echo pulse sequence will be determined by the particular 
clinical application. For prostate imaging a TE of 30 to 150 ms, a TR of 3000 to 4000 
ms and a receiver bandwidth of 16 kHz is preferred. 

[0028] Referring particularly to Figs. 5 and 6, the first step in the imaging 

method according to the present invention is to acquire two images as indicated at 
process block 400. One image is acquired with the GRE pulse sequence of Fig. 3 and 
the second is acquired with the spin-echo pulse sequence of Fig. 4. Both acquisitions 
sample the same locations in k-space to produce two 256x256 element k-space data 
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sets 401 . The two image acquisitions are done as close in time as possible so that the 
subject remains in the same position in both acquisitions. 
[0029] As indicated by process block 402, two separate images 403 are 

reconstructed from the two acquired k-space data sets 401 . This is done by performing 
a two-dimensional, complex fast Fourier transformation on each k-space data set 401 
to produce two, 256x256 pixel complex images 403. Each pixel is a complex number 
having an I and a Q value. Two corresponding phase images 405 are produced at 
process block 404 by calculating the phase ((j)) at each pixel in the two images 403: 

<j) - tan '(I/Q). 

[0030] The next step is to calculate a phase difference image as indicated at 

process block 406. The phases at corresponding pixels in each of the two phase 
images 405 are subtracted to produce a phase difference value A§ for corresponding 
pixels in the phase difference image 407: 

A<Kx,y) = A^(x,y) - A4> 2 (x,y). 

[003 1 ] Phase unwrapping may be employed at this point. A number of 

different methods well known to those skilled in the art may be used to unwrap the 
phase rolls in the phase image. A simple approach is to take orthogonal profiles 
through the center of the phase difference image resulting in two one-dimensional 
phase vectors, A(j>(x,0), and A<()(0,y) where the coordinate 0,0 represents the center of 
the image. Phase rolls are then corrected by performing the complex equation: 

e -i (A())(x,y) - A4>(x,0) - A(j)(0,y) ) _ 

[0032] As indicated by process block 408, a conventional magnitude image 

409 is also produced from the complex spin-echo image 403. Each pixel magnitude 
(M) is calculated from the corresponding complex values: 
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M(x, y) = V(I(x,y)) 2 +(Q(x,y)) 2 . 



A magnitude image can also be calculated from the complex GRE image and the 
corresponding pixel magnitudes averaged with the spin-echo magnitudes to increase 
the magnitude image SNR. 

[0033] The final step as indicated at process block 410 is to determine seed 

locations from the phase difference image 407 and superimpose their locations on the 
magnitude image 409. More specifically, the location of each brachytherapy seed is 
determined by inputting the phase difference image 407 into an automatic seed soiling 
algorithm such as that described in Brinkmann DH, Kline RW: Automated seed 
localization from CT datasets of the prostate. Medical Physics 1998; 25(9): 1667-72. 
The location of each seed is then superimposed on the magnitude image 409 since the 
phase difference and magnitude images share the same spatial coordinates. Seed 
locations within the magnitude image can then be described by modifying the 
magnitude of each pixel so that seed locations are defined by unique gray scale values, 
or by modifying the color of these pixels. Also, a graphical representation, or icon, 
with the same dimension as the true seed may be overlaid on the magnitude image at 
the proper location. 

[0034] When imaging the prostate with brachytherapy seeds implanted the 

magnitude image produced as described above differentiates between the various types 
of soft tissues. The phase difference image on the other hand clearly depicts the seeds 
as bright objects in a relatively uniform surrounding. The combined images thus 
produce an image 41 1 which reveals very clearly the bright seeds surrounded by well 
differentiated soft tissues. 

[0035] Many variations are possible from the above-described preferred 

embodiment. For example, rather than two separate pulse sequences, a single pulse 
sequence may be used which acquires both a gradient-recalled echo NMR signal and a 
spin-echo NMR signal. This is particularly useful when patient motion is a problem 
since corresponding views in both images are acquired at substantially the same time. 
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Also, while the use of a spin-echo and a gradient-echo pulse sequence is preferred, 
good results may also be achieved by acquiring two separate GRE images using 
substantially different scan parameters. For example, a pulse sequence with a TR of 
300 to 500 ms that produces two gradient recalled echoes with a TE ranging from 5 to 
10 ms for the first echo signal and from 40 to 100 ms for the second echo produces a 
good phase difference image. Also, the use of a gradient recalled echo signal pulse 
sequence to acquire two images with vastly different receiver bandwidths (e.g. 8 kHz 
and 64 kHz) will produce a good phase difference image. 

[0036] It is also possible to acquire more than one NMR signal with each pulse 

sequence. For example, multiple spin-echo signals may be acquired with a fast spin- 
echo (FSE) pulse sequence and multiple gradient-recalled echo signals may be 
acquired with an echo planar imaging (EPI) pulse sequence. It is also possible to 
acquire a single complex image using a pulse sequence in which a magnitude image 
can be reconstructed and a phase image which reveals phase perturbations at each seed 
location can be reconstructed. With current MRI systems a second phase image is 
needed as described above to distinguish such perturbations. It is contemplated that 
this may not be required in the future as MRI system phase errors are reduced. 
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